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ABSTRACT 



Aims. We present a diagnostic tool to determine the abundance of the crystalline silicate forsterite in AGB stars surrounded by a thick 
shell of silicate dust. Using six infrared spectra of high mass-loss oxygen rich AGB stars we obtain the forsterite abundance of their 
dust shells. 

Methods. We use a monte carlo radiative transfer code to calculate infrared spectra of dust enshrouded AGB stars. We vary the dust 
composition, mass-loss rate and outer radius. We focus on the strength of the 11.3 and the 33.6 fim forsterite bands, that probe the 
most recent (11.3 yum) and older (33.6 //m) mass-loss history of the star. Simple diagnostic diagrams are derived, allowing direct 
comparison to observed band strengths. 

Results. Our analysis shows that the 1 1.3 fim forsterite band is a robust indicator for the forsterite abundance of the current mass-loss 
period for AGB stars with an optically thick dust shell. The 33.6 yum band of forsterite is sensitive to changes in the density and 
the geometry of the emitting dust shell, and so a less robust indicator. Applying our method to six high mass-loss rate AGB stars 
shows that AGB stars can have forsterite abundances of 12% by mass and higher, which is more than the previously found maximum 
abundance of 5%. 

Key words. Stars: AGB and post-AGB - Stars: mass-loss - Astrochemistry - Radiative transfer - Infrared: stars 



1. Introduction 

Asymptotic Giant Branch (AGB) stars represent a late phase in 
the evolution of low- and intermediate mass stars (0.8-8.0 M ). 
At the end of the Red Giant Branch phase when the helium in 
the core is depleted a star enters the AGB phase, with a de- 
generate carbon and oxygen core surrounded by a helium and 
a hydrogen burning shell. The burning shells are surrounded 
by a convective envelope and a dynamically active atmosphere 
(Habing & Olofsson 2003). Depending on the carbon-to-oxygen 
ratio, either oxygen or carbon is locked in the CO molecule. This 
leads to two possible chemistries. Initially AGB stars exhibit an 
oxygen-rich chemistry reflecting the original composition of the 
natal molecular cloud. An AGB star becomes C-rich when suffi- 
cient carbon is dredged up from the helium-burning shell where 
it is formed by nucleosynthesis processes. In this study we will 
only consider O-rich AGB stars. 

One of the most important aspects of AGB stars is their 
stellar wind. AGB stars can have mass-loss rates between 
10 & M yr~' and about 10~ 4 M yr . The high mass-loss rate 
automatically constrains the life-time of these objects, and gives 
rise to dense circumstellar envelopes. Two mechanisms con- 
tribute to this mass-loss. First matter gets accelerated outwards 
due to shocks (Habing & Olofsson 2003). When the material 
reaches distances where the temperature is low enough the ma- 
terial can condense into solid-state particles (from here on called 
dust). Then a second acceleration mechanism occurs: radiation 
pressure on the dust grains. The gas is coupled to the dust and 
dragged along. This leads to expansion speeds of 10 km s _1 and 
higher. 



Even though amorphous silicate dust is more abundant, 
crystalline silicate dust is also seen in many astronomi- 
cal environments like disks around pre - main-sequence stars 
(IWaelkens et alj|1996t iMeeus et alJl200U iKessler-Silacci et al.l 
12006 ), comets dWoodenl 12002*1)7 post-main-sequence stars 
(IWaters et all Il996t ISvlvester et all 11999b iMolster et al l T2 002) 
and a ctive galaxies (Markwick- Kempe r et alj|2007t ISpoon et al.l 
2006). Theoretical studies predict that the formatio n of crys- 
talline silicate dust is dependent on the gas density (Tielens et al.l 
1998; iGail & Sedlmavnj 19991; ISogawa & Kozasall 19991) . Higher 
densities are more favorable for the formation of crystals than 
low densities. [Speck et al. (2008) also argue that there is a cor- 
relation between crystallinity, mass-loss rate and the initial mass 
of the star. 

These suggestions are consistent with the detection of crys- 
talline silicates in the high density winds of OH/IR stars and 
the disks around post-AGB stars. Indeed, observations show 
a lack of crystalline silicate dust bands in the spectra of low 
mass- loss rate AGB stars (I Waters et al.l 1 19961: ISvlvester et al.l 
fl999h. However, an alternative explanation for this lack exists. 
Kemper et al. (2001) have shown that a contrast effect hinders 
the detection of crystalline silicate bands in infrared spectra of 
low mass-loss rate AGB stars. This contrast effect is further dis- 
cussed in sect. 3.2.2. 



iKemper et all (120011) obtained an indication for the crys- 
talline silicate abundance in the outflow of AGB stars by com- 
paring the strength of crystalline silicate emission bands of ob- 
servations to those of model spectra. For the crystalline silicates 
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enstatite and forsterite they measured a maximum abundance of 
~5% by mass. 

AGB stars are the main contributors of crystalline silicate 
material to the ISM. The creation of crystalline silicates in the 
ISM is not possible since the temperatures are too low. Therefore 
the hypothesis would be that the abundance of crystalline silicate 
dust in the ISM would be comparable to that found in the dust 
shells of AGB stars. The crystalline silicate abundance in the dif - 
fuse ISM has been investigated by iKemper et all d2004l [2005). 
They looked at the line of sight towards the Galactic Centre and 
found an upper limit of the degree of crystallinity in the diffuse 
ISM of 2.2% by ma ss. The abundance found in AGB stars by 
Kem per et alJ (1200 ll) is therefore higher than that for the diffuse 
ISM, suggesting the crystalline material may be effectively de- 
stroyed or hidden. Possible mechanisms could be amorphization 
by inclusion of ir on atoms in the crys talline lattice, shocks and 
UV radiation (see lKemper et al.ld2004l2005l) for a discussion on 
such mechanisms). 

Determining the crystalline fraction of silicates in AGB stars, 
as a function of mass-loss rate, will enhance our understand- 
ing in two areas: First, what are the conditions needed for 
the formation of crystalline silicates? Dust condensation the- 
ory suggests that t he density is critical in the formation of 
crystalline silicates dTielens et alJ ll998T: iGail & Sedlmavrll 19991; 
Sogawa & Kozasa 1999), but this has still not been verified by 
observations. Secondly, how much crystalline silicate material is 
deposited by AGB stars to the diffuse ISM? Do the abundances 
in the outflow of AGB stars and the ISM match? 

The most thorough way to determine the forsterite abun- 
dance in the outflows of AGB stars is to fit the full spectral range 
offered by ISO or Spitzer spectra. This is an exhausting endeavor 
to apply to multiple stars and it does not necessary deliver a good 
understanding on how the spectral bands depend on dust shell 
parameter s. Therefore it is us eful to follow and expand the ap- 
proach of Ke mper et al.l d2001l) who derived diagnostics for the 
33.6 jum forsterite band. We expand on this study by carefully 
studying the dependence of the spectral bands on the mass-loss 
rate, dust composition, thermal contact between grains and we 
also use two well-chosen bands, one probing the most recent 
mass-loss (1 1.3 fim) and one probing the older, colder mass-loss 
(33.6 yum). We derive a diagnostic diagram that is fairly easy to 
use and allows one to derive the forsterite abundance in the AGB 
star outflows. 

In sect. 2 we outline how the models were computed and 
show several examples of computed spectra. The strength deter- 
mination of the two forsterite bands is discussed in sect. 3. Also 
the diagnostic diagram to measure the forsterite abundance is 
introduced and possible limitations are discussed. In sect. 4 we 
apply the diagnostic diagram to six ISO-SWS spectra of AGB 
stars. 

2. Modeling the AGB dust shell 

In this section we present the setup used to compute model spec- 
tra and show a few examples of infrared spectra. 

The model spectra were com puted using the M onte Carlo 
radiative transport code MCMax dMin et al.ll2009t) . Within the 
code we treat AGB stars as a central star enclosed in a spherically 
symmetric dust shell. The central star is defined by its radius 
and spectrum. A spectrum of an M 9III giant (2667K) is used for 
the central star dFluks et alJll994l). For the luminos ity a typical 
value of 7000 L is used dHabing & Olofssonll2003l) . Using this 
luminosity and integrating the spectrum one obtains a radius of 
395 R for the central star. 



An r~ 2 density distribution for the dust shell gives rise to a 
time-independent outflow at a constant velocity. There is ample 
evidence for strong t i me va riability of AGB mass-loss, see for 
example iDecin et all d2007l) . which would correspond to a dif- 
ferent density structure than the one adopted here. We will come 
back to this point in section 3.2.3. 

The terminal wind speed (v exp ) is set to the typical value of 
10 km s _1 . For the dust-to-gas ratio a value of 1/100 is assumed. 
For our models the dust-to-gas ratio only determines the amount 
of material in the dust shell that is in the form of dust. Changing 
this value will result in a different dust mass-loss rate and there- 
fore it only scales the total mass-loss rates used in this paper. 
The temperature of the dust shell reaches a value of 1000 K at a 
radius of ~10 AU. This radius is taken to be the radius at which 
the amorphous silicate condenses, and therefore defines the in- 
ner radius of the dust shell. Doing this we assume that the dust 
species formed before this 1000 K limit have a negligible effect 
on the spectrum. This is especially true for the optically thick 
dust shells where the infrared spectrum is not sensitive anymore 
to changes in the inner part of the dust shell. We note that in op- 
tically thin du st shells the maximu m dust temperature detected 
in a study bv lHeras & Honvl d2005l) is around 700 K. So the cal- 
culations take somewhat high condensation temperatures for op- 
tically thin dust shells. 

2.1. Dust opacities 

In this paper we only consider the major components of the dust 
in AGB outflows. We use silicate grains, both in crystalline and 
amorphous form and metallic iron grains. To calculate the opac- 
ities of dust grains, the chemical composition, lattice structure 
and shape distribution have to be chosen. These choices will be 
discussed in this section. 

The composition of the amorphous silicate component in 
AGB stars is currently not well known. In this paper we fix the 
amorphous silicate composition to be consistent with the spec- 
troscopic signature of the low mass-loss rate AGB star Mira. 
Mira shows no clear signs of crystalline silicates making it an 
ideal case to study the amorphous component in AGB outflows. 
We have fitted the ISO-SWS infrared spectrum of Mira using 
amorphous silicates with different compositions and find that the 
best fit was made with 65.7% MgFeSi0 4 and 34.3% Mg 2 Si0 4 
implying an average composition of Mgi j6Feo.64Si04. The fit 
was made with CDE shaped particles (see below). By adopting 
the composition for the amorphous silicates as found for Mira we 
make the implicit assumption that the composition of the amor- 
phous silicate dust grains is not dependent on the mass-loss rate 
of the AGB star. 

For the refractive indices we have used data measured in 
the laboratory for a large number of material compos itions. We 
use he re for the amorphous s ilicat es the data from iJager et al.l 
(I2003I) and iDorschner etail d 19951) . the forsterite data from 
Serv oin & Pirioul d!973l) . the enstatite data from Jager et al. 
(1998) and the metallic iron laboratory data from Ordal et al. 
dl988h . We always use an equal abundance of enstatite and 
forsterite. The opacities of forsterite are shown in fig.Q] 

iMin et al] d2003l) have shown that the use of homogeneous 
spheres as particle shapes introduces unrealistic effects. This is 
because the symmetries in these particles introduce resonances. 
It has been shown by severa l authors that the CD E (Continuous 
Distribution of Ellipsoids, Bohr en & Huffman! d 1983b ) shape 
distribution provides a very good fit to both astronomical obser- 
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Fig. 1. The opacities of the crystalline silicate forsterite. They 
clearly show the strong bands at 11.3 fim and 33.6 fim bands 
which are investigated in this study as possible forsterite abun- 
dance indicators. The opaciti es are calculated using th e complex 
refractive indices taken from[S ervoin &Pirioul(ll973h and using 
CDE particle shapes. 



vations and laboratory measurements. Therefore, in this study 
we use the CDE shape distribution. 

Laboratory experiments show that crystalline silicates, 
metallic iron and amorphous silicates have different con- 
densation temperatures , which are also pressure dependent 
(lGail&Sedlmavrlll999l) . It is very well possible that high tem- 
perature condensates ma y act as seeds for lower te mperature ma- 
terials to condense onto (Idail & Sedlmavd fT999>. Therefore we 
must consider chemically inhomogeneous grains as a possibil- 
ity, since this will affect their temperatures and so their IR emis- 
sion. However, grain condensation is likely to be a dynamical 
process which is difficult to model, and so it is hard to pre- 
dict from chemical equilibrium condensation theory what the 
nature of the chemically inhomogeneous grains may be. Given 
these uncertainties we decided to use two choices for chemi- 
cally inhomogeneous grains and investigate the effects of these 
two choices on the calculated spectra. These two choices should 
span most of parameter space. Our standard model assumes that 
amorphous silicate grains condense first, followed by metallic 
iron that forms on top of the amorphou s silicates, along the lines 
suggested by iGail & Sedlmavrl ([l999). Therefore metallic iron 
and amorphous silicates are in thermal contact and have the same 
temperature in our calculations. We assume the crystalline sili- 
cates to remain isolated and so not in thermal contact with other 
dust species. This minimizes the temperature of these grains, due 
to the low opacity at near-IR and optical wavelengths. Our sec- 
ond choice is to assume that all grains are in thermal contact and 
so have the same temperature. This obviously eliminates any 
temperature differences and probes the extreme case in which 
crystalline silicates, despite their low opacities at short wave- 
lengths, are still efficiently heated by stellar photons. 

2.2. Model spectra 

In this section we show examples of computed spectra. This 
will illustrate the effect of the mass-loss rate, metallic iron abun- 
dance, forsterite abundance and the outer radius of the dust shell 
on the infrared spectra. A standard set of model parameters is 




15 20 25 
Wavelength (micron) 

Fig. 2. Three model spectra with different mass-loss rates are 
shown. The solid grey, solid black and dashed black line have 
mass-loss rates of 1 • 10~ 6 M yr~', 1 • 10~ 5 M G yr _I and 1 • 
lO~ 4 M yr _1 respectively. All model parameters other than the 
mass-loss rate are set to standard values (see tableQ]). 



used (see table Q]) and for the rest of this paper, if the parame- 
ters are not mentioned explicitly, they are set to these standard 
values. 

For the metallic iron abun dance a standa r d valu e of 4% is 
chosen. This is consistent with KemperetalJ (2002), who have 
shown that the infrared spectrum of OH127. 8+0.0 can be ex- 
plained by including that amount of metallic iron in their mod- 
els. A starting value for the forsterite and enstatite abundance 
has been chosen to be 4%. In most studies the outer radius of 
the dust shell is practically set to infinity. In this study we want 
to treat the outer radius as a variable, in order to investigate its 
effect. As a standard value we have ch osen an outer radius of 
500 AU. An interesting example is that IChesneau et al I J2005h 
have shown that the dust shell of the current mass-loss of OH 
26.5+0.6 is about 250 yr old. Using a 10 km s _1 outflow veloc- 
ity this translates to an outer radius of the dust shell of roughly 
500 AU. 

Figure [2] shows three spectra with different mass-loss rates. 
These mass-loss rates are chosen to represent objects with the 
9.7 fim band in emission, partially in absorption and in absorp- 
tion. The spectrum with a mass-loss rate of 10 -6 M yr~' shows 
both the 9.7 fim and the 18 fim band of amorphous silicate in 
emission. In this case the dust shell is optically thin and the spec- 
trum of the central star can be observed. Increasing the mass-loss 
rate to 10~ 5 M yr makes the blue part of the spectrum opti- 
cally thick, obscuring the spectrum of the central star and caus- 
ing the 9.7 fim band to go into self-absorption. The spectrum of 
the 10~ 4 M yr -1 case has the 9.7 fim and the 18 fim band of 
amorphous silicate firmly in absorption. 

Spectra with different metallic iron abundances are shown in 
fig. [3] Metallic iron absorbs short wavelength radiation more ef- 
ficiently and since these grains are in thermal contact with the 
amorphous silicates, both heat up. In turn the amorphous sili- 
cates radiate this energy at longer wavelengths. It can be seen 
from the spectra that metallic iron influences the opacity in the 
3 to 8 fim region and that the emission at longer wavelengths is 
increased when metallic iron is added. 

Figure [4] shows model results for different forsterite abun- 
dances. The forsterite band at 11.3 fim can clearly be recognized 
as a shoulder on the red side of the 9.7 fim band of amorphous 
silicate, and increases in strength with increasing abundance. 

To investigate the influence of the outer radius on the spec- 
trum, models were run with outer radii up to 3000 AU (figure|5j. 
From the figure it can be seen that changing the outer radius has 
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Table 1. List of model parameters and their standard values and the range over which the parameters are varied. 



Parameter Standard value Minimum Maximum 

Mass-loss 3 • lO^Moyr" 1 1 ■ lO^Moyr 1 1 ■ l(T 4 M G yr- 

Metallic iron abundance 4% 0% 9% 

Forsterite abundance 4% 0% 10% 

Enstatite abundance set equal to forsterite 0% 10% 

Outer radius 500 AU 300 AU 3000 AU 




Wavelength (micron) 

Fig. 3. Four spectra with different metallic iron abundances are 
shown. The solid black, solid grey, dashed black and bold black 
lines have metallic iron abundances of 0%, 2%, 4% and 9% 
respectively. All model parameters other than the metallic iron 
abundance rate are set to standard values (see table[TJ). 
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Fig. 5. The effect of the outer radius of the dust shell on the spec- 
tra is plotted. The solid black, solid grey and dashed black lines 
have outer radii of 300 AU, 1000AU and 3000 AU respectively. 
All model parameters other than the outer radius are set to stan- 
dard values (see tableQ}. 
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Fig. 4. The 9.7 /im band of amorphous silicate of spectra with 
different forsterite abundances are shown. The solid grey, solid 
black and dashed black lines have forsterite abundances of 0%, 
4% and 10% respectively. All model parameters other than the 
forsterite abundance are set to standard values (see table[T]i. 

a significant impact. By increasing the outer radius there will be 
more cold dust and therefore the spectrum becomes redder. The 
outer radius also has a profound influence on the strength of the 
33.6 /mi band. For an outer radius of 3000AU the band is seen 
to be strong and easily detectable. If the radius is decreased to 
300AU the strength of the band is also seen to go down steeply. 
We will see this effect more quantitatively in sect. 3.2.2. 



We have chosen to measure the equivalent width as a measure of 
the strength of the band because it is less sensitive to a choice in 
grain shapes, than for instance the ratio of the peak flux over the 
continuum. Although the band shape may change for different 
grain shapes, the total power remains conserved. 

Both bands require a different method for measuring their 
strength. The 11.3 /mi band is an absorption band superimposed 
on the 9.7 /mi band of amorphous silicate and a careful mea- 
surement of its strength is needed. In our method to measure 
the strength of the 11.3 /mi band we limit this study to the case 
where the 9.7 /mi band of amorphous silicate is in absorption. 
As a measure of the strength of the band, the equivalent width is 
calculated from the optical depth profile (see sect. 3.1). 

The 33.6 /mi band is an emission band for the mass-loss 
rates that are considered. This means the band is predominantly 
formed in the optically thin region of the wind. Constructing a 
continuum under the band enables us to calculate the equivalent 
width. The details on how this is calculated are discussed in sect. 
3.2. 

Comparable to the 11.3 //m band of forsterite, enstatite has a 
band at 9.3 /mi. We are looking into the applicability of this en- 
statite band, but its analysis is more complicated than for the 11.3 
/mi band of forsterite. The 9.3 /mi band of enstatite lies closer to 
the absorption maximum of the 9.7 //m band of amorphous sili- 
cate. In this paper we focus on the bands of forsterite on l y, sinc e 
these can be compared with the results of Kemp eret alJ (1200 ll) . 



3. Forsterite abundance indicators 

In this section we discuss the applicability of the 11.3 /mi and 
the 33.6 /mi forsterite bands as forsterite abundance indicators. 
We compute the strength of these bands for model spectra and 
study their dependence on the model parameters. For both bands 
the strength is measured in spectra with varying mass-loss rate, 
forsterite abundance, metallic iron abundance and outer radius. 



3.1. The 11.3 //m forsterite band 
3.1.1. Method 

To measure the strength of the 11.3 //m band a local continuum 
is defined on top of the 9.7 /mi amorphous silicate band. A linear 
interpolation between the flux densities at 8.0 /zm to 8. 1 /zm and 
12.9 //m to 13.0 /mi is made to construct this local continuum. 



B.L. de Vries et al.: Determining the forsterite abundance in AGB stars 



5 





1 10" 


M e yr 1 




/ 5 10" 


M e yr _1 "%\ 




/ 2 10" 











10 11 
Wavelength (micron) 




11.0 11.2 11.4 
Wavelength (micron) 



Fig. 6. Three optical depth curves of the 9.7 fim band of amor- 
phous silicate for three different mass-loss rates. The constructed 
local "continuum" in order to extract the 11.3 fim band is also 
shown. The curves are spaced vertically for clarity. The bottom, 
middle and top curve have mass-loss rates of 2 • 10~ 5 Moyr _1 , 
5 ■ 10~ 5 M G yr~' and 1 ■ 10~ 4 M G yr~' , respectively. All model pa- 
rameters not mentioned in the graph are set to standard values 
(see table [TJ 



For every wavelength in the range of the 9.7 /mi band the 
optical depth is calculated: F(A)/F conl (A) = e~ T , where F (A) is 
the flux of the model and F co „,(A) the flux of the continuum. 
In the optical depth profile t(A) a local "continuum" is defined 
under the 11.3 jum band. This is done by making a linear fit 
through the optical depths at the wavelength points 10.5 fim to 
10.6 fim and 11.6 fim to 11.7 fim. Figure [6] shows three opti- 
cal depth curves together with the continuum constructed under 
the 11.3 fim band. The 11.3 fim band can now be normalized: 
T(A)/T cont (A) - 1. The equivalent width in terms of optical depth 
is calculated by numerical integration from 10.5 fim to 1 1.7 fim. 
Figure [7] shows three normalized optical depth curves for the 
11.3 fim band. 

3.1.2. Results 

The effect of changes in the mass-loss rate on the strength of 
the 11.3 fim band is shown in fig. [8] for a range of forsterite 
abundances. The band strength decreases with increasing mass- 
loss, reaching a constant value for the highest mass-loss rates. 
This behavior can be understood as follows. For the highest 
mass-loss rates, optical depths are high for many lines of sight 
through the dust shell. In this case the band strength is deter- 
mined by the ratios of the optical depths of the dust species, 
and is not sensitive any more to temperature differences. This is 
why the curves in fig. [8]reach a constant value. At lower mass- 
loss rates, emission from optically thin lines of sight through the 
dust shell become more important. The temperature difference 
between dust species plays a role and the warmer amorphous 
silicates and metallic iron will begin to fill in the 9.7 fim ab- 
sorption band seen so prominently at high mass-loss rates. This 
results in a relatively smaller optical depth of the amorphous sil- 
icate band compared to that of the crystalline silicates, and so the 
ratio (T(A)/T cont (A) - 1) will increase. The curve for 2% forsterite 
does not follow the general trend seen in fig. [8] at low mass-loss 
rates. This is due to the approximate way in which the contin- 



Fig. 7. Three normalized optical depth curves of the 11.3 fim 
band of forsterite for three different forsterite abundances. From 
the solid, dashed, dash-dotted curves have a forsterite abundance 
of 2%, 4% and 10% respectively. All model parameters not men- 
tioned in the graph are set to standard values (see table [TJ 
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Fig. 8. The equivalent width of the 11.3 fim band of forsterite as a 
function of the mass-loss rate. The different curves have different 
forsterite abundances. All model parameters not mentioned in 
the graph are set to standard values (see table [TJ. The upward 
trend seen for the 2% forsterite curve in the mass-loss rate range 
of 1 ■ 10~ 5 M yr~' to 5 • 10~ 5 M yr -1 is due to the approximate 
way in which the continuum is estimated in order to calculate 
the strength of the 11.3 fim (see section 3.1.2). 



uum is estimated (linear fit), causing small errors that can be 
important for weak bands. 

The outer radius of the dust shell has no significant effect on 
the strength of the 11.3 fim band and a figure has been omitted. 

Adding metallic iron increases the emission of the contin- 
uum under the 9.7 fim band. This reduces the optical depth in 
the region of the 9.7 fim band (T conr (A)), causing the equivalent 
width of the 1 1.3 fim band to increase (see fig.[9j. 

In order to test the effects of our choice regarding thermal 
contact we also consider the case of thermal contact between all 
grains (section 2.1). FigureQJJcompares the strength of the 1 1.3 
fim band in the case where only amorphous silicate and metallic 
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Fig. 9. The effect of changes in the metallic iron abundance on 
the strength of the 11.3 fim forsterite band for two choices of 
the forsterite abundance of 2% and 10% respectively. Solid lines 
represent 2% metallic iron, dashed lines 10%. All model param- 
eters not mentioned in the graph are set to standard values (see 
table[Qi. 
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Fig. 10. The equivalent width of the 11.3 fim band is plotted as 
a function of the mass-loss rate. It shows the 2% and the 10% 
forsterite abundance curves for two different cases. The first case 
is where only metallic iron and amorphous silicate grains are in 
thermal contact (solid line). In the second case all dust species 
are in thermal contact (dashed line). All model parameters not 
mentioned in the graph are set to standard values (see table[TJ. 

iron are in thermal contact with the case where all dust species 
are in thermal contact. The strength of the 11.3 fim band de- 
creases by adopting thermal contact between all the dust species. 
This is because part of the thermal energy is redistributed to the 
crystalline silicate species, which do not contribute to the con- 
tinuum. This causes the flux of the continuum to drop, resulting 
in a weaker 11.3 fim band. 

3.1.3. Discussion 

The 11.3 fim forsterite band is formed close to the star and 
so it probes the most recent mass-loss rate. We find that its 



Fig. 11. The continuum constructed under the 33.6 fim band of 
forsterite for a model spectrum with standard parameters (see 
table [D 



strength is relatively insensitive to model parameters (other than 
the forsterite abundance itself) when the mass-loss rate exceeds 
roughly 3 ■ 10~ 5 M yr~'. The largest uncertainty at these high 
mass-loss rates is due to changes in the metallic iron abundance. 
The effect of thermal contact between forsterite and other dust 
species is minor (1-2%). In our standard model we assume no 
thermal contact, which leads to conservative estimates of the 
forsterite abundance. 

3.2. The 33.6 fim forsterite band 

3.2.1. Method 

In order to measure the equivalent width of the 33.6 fim band 
a local continuum is defined by making a linear fit through the 
flux densities at 31.0 fim to 31.5 fim and 35 fim to 35.5 fim. An 
example of such a continuum is shown in fig. QT| This contin- 
uum is used to normalize the spectrum in this region, calculating 
F(X)I F cont (X) - 1. Three such normalized bands with different 
forsterite abundances are shown in fig. [12] The equivalent width 
of the normalized band is calculated by numerical integration 
between the wavelength values of 31 fim and 35.5 fim. We note 
that due to the emission nature of the band its equivalent width 
will be sensitive to the abundance and temperature of the dust 
species emitting at 33.6 fim. 

3.2.2. Results 

Figure Q~3] shows the equivalent width of the 33.6 fim band 
as function of mass-loss rate for different forsterite abun- 
dances. The mass-loss rates range from lO _6 M yr _1 up to 
2 • lO _4 M yr _1 . The strength of the band initially increases 
with increasing mass-loss rate, followed by a drop and eventu- 
ally the band is in ab s orptio n. This behavior is similar to that 
seen by Kem per et alJ yOOl) and has been explained by these 
authors in terms of changes in the temperature of the different 
dust species as a function of mass-loss rate. For low mass-loss 
rates, the dust shell is optically thin and the (near-IR) stellar pho- 
tons heat the dust directly. This causes large temperature differ- 
ences between amorphous silicate (which has high opacity at 
near-IR wavelengths) and crystalline silicates (which are almost 
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Fig. 12. Normalized 33.6 fim bands for different forsterite (fo.) 
abundances. Other parameters are set to standard values (see ta- 
ble[B 
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Fig. 13. The equivalent width of the 33.6 /mi band of forsterite as 
a function of the mass-loss rate. The different curves have differ- 
ent forsterite abundances. All model parameters not mentioned 
in the graph are set to standard values (see table[T]l. 



transparent in the near-IR). The band strength of forsterite will 
thus be weak. For higher mass-loss rates the dust shell will be- 
come more optically thick and more and more of the dust will 
receive longer wavelength photons emitted by the optically thick 
inner dust shell. The opacity of amorphous and crystalline dust 
at mid-IR wavelengths is very comparable and so the temper- 
ature difference between them will decrease. This results in an 
increase of the band strength. For even higher mass loss rates 
(larger than about 1CT 5 M Q yr~') some lines of sight become op- 
tically thick at 33.6 fim and this will result in a weakening of the 
band, eventually even leading to absorption. 

The fact that the curve for 2% forsterite becomes negative at 
low mass-loss rates is due to the construction of the continuum, 
which is curved but approximated by a linear fit. This introduces 
small errors that can be important for weak bands, especially at 
low mass-loss rates. 

The effect of metallic iron on the equivalent width of the 33.6 
yum band is demonstrated in fig. [14] Metallic iron increases the 
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Fig. 14. The equivalent width of the 33.6 //m band as a function 
of the mass-loss. Both the curves with 2% and 10% forsterite 
are shown for the metallic iron abundances of 2% (dashed line) 
and 10% (solid line). All model parameters not mentioned in the 
graph are set to standard values (see table [TJ. 
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Fig. 15. The equivalent width of the 33.6 fim band as a func- 
tion of the mass-loss rate. It shows the 2% and the 10% 
forsterite abundance curves with outer radii of 500 (solid line) 
and 2000AU (dashed line). All model parameters not mentioned 
in the graph are set to standard values (see table[TJ. 



flux of the continuum, but not that of forsterite. Since the equiv- 
alent width is a measure of the relative strength of forsterite to 
the continuum, increasing the metallic iron abundance will de- 
crease the equivalent width. The reason for this is that forsterite 
is considered not to be in thermal contact with the amorphous 
silicates, but metallic iron is. 

Increasing the outer radius of the dust shell increases the 
amount of cold, emitting dust and therefore increases the emis- 
sion of the 33.6 fim band (figure [TBIl. It is interesting that the 
effect of the size of the dust shell becomes smaller and even neg- 
ligible for low mass-loss rates. Below ~ 6 • lO~ 6 M yr~' the 
effect is already very small. The reason is that the wind is so 
optically thin that adding more material by increasing the outer 
radius adds material of such low temperatures that its emission 
is negligible. 
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Fig. 16. The equivalent width of the 33.6 fim band is plotted as 
a function of the mass-loss rate. It shows the 2% and the 10% 
forsterite abundance curves for two different cases. The first case 
(solid line) represents models in which only the metallic iron and 
amorphous silicate grains are in thermal contact. In the second 
case (dashed line) all the dust species are in thermal contact with 
each other. All model parameters not mentioned in the graph are 
set to standard values (see tableQ}. 



Figure [T6l shows the equivalent width of the 33.6 /mi band 
for the two different cases of thermal contact (see sect. 3.1.2) 
for models with 2% and 10% forsterite. It shows that the equiv- 
alent width of the band goes up when all the dust species are in 
thermal contact. In the case where all dust species are in thermal 
contact with each other, forsterite is heated up by metallic iron, 
creating a stronger emission band. The strength of this effect de- 
pends on the forsterite abundance. 



3.2.3. Discussion 

The 33.6 fim band is sensitive to the mass-loss rate of the AGB 
star as well the outer radius of the dust shell. The size of the dust 
shell has an effect because the 33.6 fim band is mainly formed by 
the emission of cold dust in the outskirts of the wind. This tells us 
that other parameters that influence this cold dust will also affect 
the 33.6 /mi band. As an example we know that AGB stars can 
have periods with different mass-loss rates and material of such 
a previous mass-loss period could also influence the 33.6 /mi 
band. Such time variable mass-loss can be modeled as a sudden 
jump in the density, but can also be studied by introducing a 
steeper or flatter radial density gradient. Changes in the radial 
density gradient will affect the relative strength of the 33.6 /mi 
emission band with respect to that of the 11.3 /mi absorption 
band. We note however that the 1 1.3 fim band is most sensitive 
to the innermost regions of the shell since it is in absorption. 
Therefore the 11.3 /mi band remains a robust indicator of the 
most recent mass-loss rate. But because of the dependence on the 
dust shell parameters the 33.6 fim band is not a good forsterite 
abundance indicator. 

Another reason why the 33.6 fim band is not a practical indi- 
cator is that for very high mass-loss rates, typically bigger than 
10~ 4 M Q /yr, all model curves converge and so become insensi- 
tive to the forsterite abundance. 



Table 2. A selection of AGB stars with their mass-loss rates 
and the results. It lists the strength of the 11.3 fim band to- 
gether with the indicated forsterite (fo.) abundance using the 
11.3 f im band. The sour ces have been selected using the cata- 
log of Sloan et al. (2003). The forsterite abundance obtained by 
iKemper et al.l (12001 ) has also been included. 



Source 


M 

(Mo yr 

(a ) 


>) 


Strength 
1 1.3 fim 
band 


Fo. abun. 
1 1.3 fim 
band 


Fo abun. 
Kemper et 
al. (2001) 


IRAS 17010-3840 


? 




0.14 


>12% 




OH104.91+2.41 


5.6 10- ; 


(1) 


0.03 


<2% 


<5% 


OH26.5+0.6 


16 10~ 5 


(2) 


0.06 


-4% 


<5% 


OH127.8+0.0 


20 10" 5 


(3) 


0.04 


-2% 


<5% 


OH30. 1-0.7 


15 10 s 


(4) 


0.10 


-9% 




OH32.8-0.3 


16 10' 5 


(5) 


0.14 


>12% 


<5% 



1: iRiechers et alJ |2005h . 
[ Justtanont & Tielensl jl992IK 4: 
iGroen ewegen ( 1 9940 



2: lJusttanont et ail 
ISchutte & Tielensl 



4. Measuring the forsterite abundance 

We apply our method to six spectra of hi gh mass-loss rate AGB 
stars (see table© observed by ISO-SW S (ISvlvester et"aT][T 999). 
The reduced data have been taken from lSloan et al.l d2003l) . We 
have used ISO-SWS spectra and not Spitzer spectra, since those 
are still in the pipe-line to be reduced and analyzed. Inspection of 
the spectra shows that the 9.7 fim amorphous silicate band is in 
absorption. In addition, a shoulder is seen in the observed spectra 
at 11.3 fim, which we attribute to forsterite. The lates t detection 
of a ve ry strong 1 1.3 fim band has been reported by Spec k et aD 
(T2Q08h . Their detection is the first detection where the 11.3 fim 
forsterite band does not appear as a shoulder on the 9.7 fim band 
of amorphous silicate, but as a separate absorption band. Given 
the presence of the 33.6 fim forsterite band in the ISO spectra of 
our sample of six stars dKemper et al.l l2001). we suggest that the 
1 1.3 fim shoulder can indeed be attributed to forsterite as well. 

We have also fitted the 9.7 fim band of amorphous silicate 
using the opacities discussed in section 2. 1 and find good agree- 
ment (not shown). This suggests that the adopted opacities that 
were obtained by fitting the 9.7 /mi band of Mira also apply to 
the higher mass-loss rate AGB stars. It also means that we can 
use these opacities in this section to measure the forsterite abun- 
dance. 

We have measured the band strength of the 11.3 fim forsterite 
band as outlined in sect. 3.1. The errors calculated for the 
strength of the bands are determined using standard error prop- 
agation formula ap plied to the errors in the flux obtained from 
Sloan et al] (120031) . A 10% error in the flux of the continuum is 
introduced to account for possible errors in the construction of 
the continuum. This is done for the two continua constructed in 
order to measure the strength of the 11.3 fim band. Mass-loss 
rates from the literature are listed in Table [2] No mass-loss rate 
was published for IRAS 17010-3840. The spectrum of IRAS 
17010-3840 has its 9.7 fim band of amorphous silicate firmly in 
absorption (similar to the other sources), therefore its mass-loss 
rate has been assumed to lie between 5 • 10 -5 M yr -1 and 2 • 10~ 4 
M Q yr 1 . 

Figure [T71 shows the strength of the 11.3 fim band for the 
six sources. The forsterite abundances indicated by fig. [T7] are 
listed in table [2] The results show that both stars with a high 
and low forsterite abundance are found (below 2% and up to 
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Fig. 17. The strength of the 11.3 band for six ISO-SWS 
sources (see table __} together with the model results from fig. 
[8] All model parameters not mentioned in the graph are set to 
standard values (see table [1_> 

~12%). Because of our choice of thermal contact between the 
dust species these values are lower limits. 

Kem per et alJ d200ll) have used the 33.6 fjm band to mea- 
sure the forsterite abundance in AGB stars. From their sam- 
ple of six AGB stars (including OH104.91+2.41, OH127.8+0.0, 
OH26.5+0.6, OH32.8-0.3) they found a maximum forsterite 
abundance of 5% by mass. The forsterite abundances we find 
for OH104.91+2.41, OH127.8+0.0 and OH26.5+0.6 do not ex- 
ceed this maximum of 5%. The sources OH32.8-0.3, OH30.1- 
0.7 and IRAS 17010-3840 have significantly higher forsterite 
abundances. 

5. Conclusions 

We have shown that for optically thick dust shells the strength 
of the 11.3 fj.m band is a robust quantity to use as an indicator of 
the forsterite abundance of the current mass-loss period. 

Applying the 33.6 //m band is practically impossible without 
knowledge of the mass-loss rate, the size of the dust shell and 
other parameters that influence the emission of the cold dust in 
the outskirts of the dust shell of the AGB star. Additionally, for 
certain mass-loss rates (depending on the dust shell parameters) 
the 33.6 fim band is also in transition from an emission to an 
absorption band, making the use of this band impossible. 

As an example and first result we used the 11.3 //m band as 
an indicator of the forsterite abundance for six AGB stars. This 
showed that these objects can have low (below 2%) but also very 
high forsterite abundances (12% and possibly higher), mean- 
ing that AGB stars can have a forsterite abundan ce higher than 
5% by mass, which was the maximum found by IKemper et alJ 
J200ll) . It also means that the discrepancy between the amount 
of crystalline material in the ISM and that in AGB stars could 
even be higher. 

If the forsterite abundances found for the six sources in this 
study represent all AGB stars, the forsterite abundance in AGB 
stars alone would be higher t han the crysta l line a bundance (up- 
per limit of 2.2%) found bv IKemper etaT] (120051) for the ISM. 
To find out if the abundances found for these stars are typical 
for AGB stars, more spectra have to be analyzed. The number 
of spectra are also not enough to test any correlations. But this 



paper opens the way to a more elaborate study similar to the one 
described in sect. 4 of this paper. Determining the amount of 
forsterite that is produced by high mass-loss rate AGB stars will 
help to test the suggested correlations and see if there is really a 
discrepancy between the produced abundance in AGB stars and 
the abundance found in the ISM. 
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